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The  phonon  relaxation  of  the  vibrational  adbond  of  an  adsorbed  no lecule 
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and  a  phonon -daaped  adbond  Irradiated  by  a  laser  are  studied.  In  the  first 

half  of  the  paper,  approx laat ions  are  Bade  within  the  Zwanzlg  projection 

operator  forasllsa  in  order  to  arrive  at  a  aaster  equation  for  the  reduced 

density  operator  of  a  Basil  sub ays tea  (the  adbond)  in  contact  with  a 

reservoir  (the  phonons).  The  conditions  of  validity  for  the  Born  and  Markov 

approx laat ions  are  derived.  It  is  shown  that  the  aaster  equation  is  only 

valid  for  tlaes  t  »  t  ,  where  t  Is  the  characteristic  tias  of  the 

c  o 

reservoir.  These  results  are  then  applied  to  the  phonon  relaxation  of  the 
vibrational  adbond  of  physlsorbed  aolecules.  It  is  shown  that  for  CO 
adsorbed  on  Mi  or  Cu  (i  strongly  bound  physlsorbed  systea)  the  Born  and 
Markov  approx  last  ions  are  not  Justified.  For  the  weakly-bound  systea  A r  on 
V,  nuaerical  results  show  that  these  approx  laat  ions  can  be  aade.  Finally, 
an  adbond  interacting  with  both  laser  radiation  and  lattice  vibrations  is 
oonsidered.  This  systea  can  be  regarded  as  a  subsystea  (the  adbond)  In 
oontact  with  two  reservoirs,  where  the  oondltlons  for  validity  of  the  Markov 
approxiaation  Is  then  seen  to  be  aore  severe  than  when  each  reservoir  is 
oonsidered  independently.  For  the  phonons,  these  conditions  can  never  be 
aatched.  However,  for  an  initial  state  given  by  an  adbond  in  equlllbrlua 
with  the  lattice  vibrations,  the  oondltions  for  validity  of  the 
approx laat ions  prove  to  be  the  saaa  as  for  the  phonons  and  the  laser 
oonsidered  independently. 
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1.  Introduction 

In  quantua  mechanics,  complete  information  about  a  syataa  is  contained 
In  the  wave  function  |+>.  When  there  la  only  a  Halted  knowledge  about  a 
ayatea,  the  proper  way  to  deacrlbe  lta  atate  la  by  the  denalty  operator, 

P  •  Ip* (♦><♦(  (D 

♦ 

where  the  real  nuabera  p^  fora  a  aet  of  probability  aaplltudes. 

Ip*  -  1.  (2) 

♦ 

All  inforaatlon  about  the  propertiea  of  a  subsystem  is  given  by  the  reduced 
denalty  operator, 

•  •  Tr  p,  (3) 

where  Tr  denotes  the  trace  over  all  variables  other  than  those  of  the 

aubaystea. 

from  the  Sehrddlnger  equation,  an  equation  for  the  tlae  evolution  of 
the  density  operator  of  a  closed  system  Is  easily  derived  to  be 

1M*dt^  *  <*> 

tdiere  H  la  the  Haailtonlan  of  the  system.  From  Eq.  (A)  an  expression  for 
the  tlae  derivative  of  the  reduced  density  operator  is  obtained  by  taking 
the  trace  on  both  aides.  However,  this  tlae  derivative  still  depends  on  the 
eoaplete  denalty  operator,  becauae  the  operators  Tr  and  H  do  not  eoaaute. 

It  follows  that  to  obtain  the  tlae  evolution  of  e(t)  exactly,  the  equation 
of  aotlon  for  the  eoaplete  denalty  operator  [Eq.  (A)]  need  to  be  solved. 

Often,  we  are  interested  In  the  properties  of  a  very  aaall  ayatea 
(e.g.,  one  aolecule)  in  oontact  with  a  much  larger  one,  called  the 
reservoir.  We  assuae  that  the  reservoir  can  be  aade  arbitrarily  large,  such 


that  there  is  no  appraclabla  dlffaranca  batwaan  the  reduced  density  operator 
for  the  reservoir,  Tr#p,  and  the  density  operator  for  the  freely  in  tlae 
evolving  reservoir,  p®.  Zt  is  assuaed  that  the  latter  is  known.  When  the 

D 

interaction  is  weak,  we  can  try  to  obtain  an  approximate  solution  for  the 
tine  evolution  of  the  reduced  density  operator  of  the  subsystem.  If  the 
interaction  is  taken  into  account  to  lowest  order  in  perturbation  theory,  it 
is  easily  shown  that  a  first-order  differential  equation  is  obtained  for  the 
reduced  density  operator.1  Apart  from  o(t),  it  involves  only  the  zeroth- 
order  approximation  to  the  reservoir  density  operator,  p®.  For  longer 
tines,  the  effect  of  the  interaction  tends  to  build  up.  Therefore,  there 
will  be  a  tine  beyond  which  perturbation  theory  is  not  valid.  To  obtain  an 
equation  of  notion  for  the  subsysten  reduced  density  operator  for  all  tines, 
a  second  approx inat ion  is  necessary,  usually  denoted  as  the  Harkov 
approximation.1 

One  nethod  for  obtaining  an  equation  of  notion  for  the  subsysten 

2  3 

reduced  density  operator  is  the  Zwanzig  projection  technique.  ’  With  this 

technique,  an  exact  lntegro-dlfferentlal  equation  for  the  reduced  density 

operator  of  the  subsysten  is  obtained.  The  two  approxiaations  necessary  to 

arrive  at  a  slaple  first-order  differential  equation  for  o(t)  are  conmonly 

denoted  as  the  Born  approximation  and  the  Harkov  approxlaatlon.*  Recently, 

this  nethod  has  been  used  to  obtain  the  effect  of  the  lattice  vibrations  of 

a  orystal  on  the  dynaaics  of  an  adsorbed  aton.*  Zn  the  derivation  of  the 

■aster  equation,  the  Born  approximation  and  the  Harkov  approximation  have 

been  introduced  as  separate  approxiaations.*’*  In  later  studies  the 

6-1 1 

validity  of  both  approxiaations  have  been  studied  independently. 

Zn  this  paper  we  shall  consider  the  validity  of  both  approxiaations 
flron  a  aore  fundamental  point  of  view.  Zn  Section  2  we  give  a  review  of  the 


Z  vans  If  projection  operator  formalism.  In  the  next  two  sections  the  Born 

end  Harkov  approx  last  Ions  will  he  Introduced,  and  the  conditions  will  be 

derived  under  which  they  are  Justified.  It  will  be  shown  that  the  condition 

for  validity  is  the  aaae  for  both  approx  last  ions.  Another  widely-used 

technique  for  obtaining  an  equation  of  notion  for  the  reduced  density 

1  12 

operator  of  the  subsystea  is  the  so-called  reservoir  theory.  ’  In  Section 
5  the  Zwanzig  projection  technique  will  be  coapared  with  reservoir  theory. 

In  particular,  the  correspondence  of  the  two  approximations  involved  will  be 
discussed.  In  Sections  6-8  we  shall  discuss  the  systee  of  a  vlbratlonally 
daaped  adatoa  and  of  an  adatoa  irradiated  by  laser  light  in  wore  detail. 
Numerical  examples  will  be  given  based  on  previous  experimental  and  theore¬ 
tical  work.  Our  conclusions  are  presented  in  Section  9. 

2.  Zwanzig  projection  operator  technique 

We  consider  a  small  subsystem  In  contact  with  a  large  reservoir.  He 

denote  the  Hamiltonian  of  the  subsystea  by  H  ,  of  the  reservoir  by  H.  ,  and 

S  D 

of  their  Interaction  by  V,  such  that  the  Hamiltonian  of  the  complete  system 
is  given  by, 

H  •  Hs  ♦  «b  ♦  v-  (5 ) 

He  then  define  the  projection  operator  P  by, 

*  *  Pb  ^b*  (6> 

where  is  any  density  operator  of  the  reservoir  and  Tir^  denotes  the  trace 
over  the  reservoir  variables.  In  the  Zwanzig  projection  technique,  an 
equation  of  motion  for  Pp(t)  is  derived,  from  which  an  equation  of  motion 
for  the  subsyatem  reduced  density  operator  is  directly  obtained  by  taking 
Tr^.  The  starting  point  is  the  equation  of  motion  for  p(t),  Eq.  (*).  Let  L 


b«  the  Llouvllle  operator  defined  by, 

IP  -  gCH.p],  (7) 

with  analogous  definitions  for  L#,  L^,  Lj  and  Lfi  •  L#  ♦  L&.  Then,  fron  Eq. 
(*>,  an  exact  equation  for  Pp  can  be  derived  as^ 

♦  PL  expC -1( 1-P)Lt3(1-P)p(0) 
t 

-1  j dt *  PL  expC -i( 1-P)Lt*](l-P)LPp(t-t* ).  (8) 

0 

The  first  tern  in  Eq.  (8)  has  a  contribution  of  zeroth-  and  first-order  in 
the  interaction,  and  the  last  tern  has  contributions  of  all  orders  in  V.  To 
separate  the  free  evolution  of  the  subsysten  (deternined  by  L  )  fron  the 
interaction  with  the  reservoir,  the  following  conditions  are  lnposed  on  p& 
and  V:^ 


cvV  ’  [vV  •  °- 

(9) 

'"b  •  °- 

(10) 

Equation  (9)  puts  restrictions  on  p^,  and  Eq.  (10)  defines  V  and  H#  for  a 
given  pfc .  The  equation  of  notion  for  Pp(t)  now  becones  [Eq.  (8)] 

•  •jjj[H#,Pp(t)]  ♦  jPLIexp[-i(1-P)Lt](1-P)p(0) 
t 

-  jdt'  PLjexp[-i( 1-P)Lt'3(1-P)LjPp(t-t* ).  (11) 

0 

The  first  tern  in  Eq.  (11)  describes  the  free  evolution  of  Pp(t),  whereas 
the  last  tern  only  gives  corrections  of  second-  and  higher-order  in  the 
interaction. 

To  be  able  to  estlnate  the  contributions  of  the  different  terns,  we 
assune  that  V  can  be  written  as  a  sun  of  products  of  a  subsysten  operator 


i  yp 

and  a  rasarvolr  operator.  *  for  only  on*  tarn  we  have 

V  •  *  S.  (12) 

Let  us  define  the  operator  R(t)  in  the  interaction  picture  by1 

(*(t)  •  expert]  R  exp[-jHbt].  03) 

Furthermore,  we  assuae  that  the  reservoir  has  a  characteristic  tine  t  , 

c 

given  by 

TrblR(0)R(t)p°}  -  Trb{R2(0)p®)exp:-|t|/te3.  O*) 

It  can  then  be  shown  that  the  contribution  of  p(0)  (second  term  in  Cq.  (11)) 

to  Pp(t)  Is  a  factor  i  /t  smaller  than  the  contribution  of  the  third  tern, 

c 

If  Ap(0)  ■  p (0 )  -  pb«(0)  Is  small.  Therefore,  the  second  term  can  be 

neglected  if  we  are  only  interested  in  the  time  evolution  of  the  subsystem 

for  times  t»t  . 

c 

If  Ap(0)  is  not  small,  the  contributions  to  Pp(t)  can  be  much  larger 
than  quoted  above.  Requiring  Ap(0)  to  be  small  leads  to  another  condition 
on  pb: 

pb  *  TrV(t)‘  05) 

For  a  large  reservoir  In  thermal  equilibrium,  we  can  use 

»b  •  s  ■  —  -ir-  (,‘> 

V  b 

Ve  note  that  pb  -  Trgp(0)  is  not  a  good  choice  because  of  Eq.  (9). 

3.  horn  approximation 

The  Born  approximation  Involves  the  retention  of  only  the  terms  of 

4 

lowest  order  In  V  in  Eq.  (11).  These  are  obtained  by  replacing  L  by  LQ  in 
the  exponential  function.  The  resulting  equation  for  the  subsystem  reduced 


density  operator  Is 


t 

-  *nf£H8,o(t)]  -  Trjdt*  LjexpC-lLQfDLje^U-t*).  (17) 

0 


In  order  that  this  be  a  good  approx  last  Ion,  the  tern  of  next  order  In  V  oust 
be  such  sea ller  than  the  tern  of  second  order.  To  test  the  validity  of  this 
approximation  we  need  an  estimate  for  the  value  of  both  terms. 

We  assume  that,  following  Eq.  (10),  Trbv^Pb  •  0,  then  the  lowest-order 
correction  to  the  Born  approximation  is  of  fourth  order  in  V.  That  is,  we 
need  the  second-order  terms  In  the  expansion  of  the  exponential  function. 
Neglecting  the  commutator  of  Lj  with  Lq,  these  second  order  terms  are  given 
by 

♦(-if)2(1-P)LIexp[-i(1-P)L0t'](1-P)L1.  (18) 

do 

The  fourth-order  contribution  to  —  becomes 

at 

t 

•  Trbj  dt’  Li  4  expt-lL0t*]Lj  LjPbo(t-f).  (19) 

0 


Use  of  Eqs.  (12)  -  (1*0  shows  that  the  integrand  involves  a  reservoir 
correlation  function  Trb(R2(0)R2(-t» )} ,  which  is  approximately 

Trb{PbR2(°)R2(-t,)}  *  CTrb(pbR2)]2exp[-|t’|/Tc].  (20) 


The  order  of  magnitude  of  the  subsystem  factor  in  the  Integrand  is  given  by 

(2D 


lw  t '  lw  t  ' 

Tr^sVoMe  0  *  CTh#(S2o(0)]2e  0  , 


where  la  a  typical  transition  frequency  of  the  subsystem.  Given  the 


dief  Ini  t  ion 
2 


V-  -  Trb{pbR2)Tr#{o(0)S2), 


(22) 


the  Integral  of  Eq.  (19)  can  be  evaluated  by  using  these  approximations  to 


yield 


^do (t) j(4 ) 
dt 


v(t>* 


(23) 


Hereby  is  the  upperbound  of  the  integral  extended  from  t  to  infinity,  which 

Te 

because  of  Eq.  (20)  Introduces  an  error  of  the  order  (— ),  which  la 

comparable  to  the  error  introduced  by  neglecting  the  contribution  of  p(0). 

le  f 

The  effect  of  the  factor  e  In  Eq.  (21 )  la  to  alngle  out  the  Fourier 

teraa  in  the  exp&nalon  of  the  correlation  function  with  u2^.  An  analogous 

1 2 

calculation  for  the  second-order  tern  gives 


cdq(t>3(2) 

dt 


Pbo(t). 


(24) 


For  the  Born  approxleation  to  be  valid  we  have  the  condition 
^(A)  k<  ^(2)^ 


(25) 


where  Y^n*  la  the  n-th-order  contribution  In  the  interaction  to  the 
relaxation  constant,  Y.  The  latter  is  defined  as  a  typical  eatrlx  element 
of  the  Llouvllle  operator  describing  the  time  evolution  of  o(t)  due  to  the 
interaction  with  the  reservoir  (second  term  in  Eq.  (17)).  Substitution  of 
Eqs.  (23)  and  (24)  into  (25)  gives 


«  1. 


(26) 


b.  Markov  approximation 

Equation  (11)  Mows  that  the  time-derivative  of  o(t)  depends  on  the 
▼alue  of  o(t)  on  previous  times.  In  the  Markov  approxleation  this  lntegro- 
dlfferentlal  equation  la  replaced  by  a  simple  first-order  differential 
equation  for  e(t).  From  Eq.  (11)  it  follows  that  o(t)  is  a  fast -oscillating 


function,  due  to  the  free  evolution  of  the  subsystem.  These  rapid 
oscillations  can  be  ellalnated  by  transformation  to  the  interaction 
picture,1 

B(t)  -  exptj-H#t]  o(t)  exp[-jjHat].  (27) 

The  equation  of  motion  for  3(t)  Is  [Eq.  (11)] 

t 

-  -expC^t  ]TrJdt'  {Ljexpt-K  1-P)Lt* J(l-P)LIP5exp[-±Hs(t-t’ )] 

0 

b(t-f)exp[^Hs(t-t']}exp[-iHat].  (28) 

In  the  previous  section  it  was  derived  that  the  integrand  involves  a  factor 
exp[-| t» | /tc] .  Therefore,  if  3(t)  does  not  change  much  on  a  time  scale  of 
the  order  t  ,  8(t-t*)  may  be  replaced  by  3(t)  in  the  integrand.  The 
temporal  changes  in  3(t)  are  of  the  order  Y.At,  so  that  we  obtain  as  the 
condition  for  the  validity  of  the  Markov  approximation 

y  «  t"1.  (29) 

c 

y  Is  already  approximately  evaluated  in  the  previous  section  [Eq.  (2*0]. 

The  condition  under  which  the  Markov  approximation  is  Justified  then  becomes 

•\2 

— ~  «  1.  (30) 

tr 

Comparison  of  Eq.  (26)  and  Eq.  (30)  shows  that  the  conditions  for  validity 
of  the  Born  approximation  and  of  the  Markov  approximation  are  the  same.  A 
similar  result  was  obtained  by  Fano1^  In  his  treatment  of  pressure 
broadening. 

Within  the  Born  and  Markov  approximations,  the  equation  of  motion  for 
the  reduced  density  operator  of  the  subsystem  is 

■  TRcV,(t,J  -  r,(t)- 


(3D 


where  r  la  tha  Llouvllle  operator  describing  relaxation 

T  -  Trjdt’  LjexpC  — li-Qt * jLjp^expC  1  Lgt ' ] .  (32) 

0 

Substitution  of  the  solution  of  Eq.  (3D  into  Eq.  (ID  shows  that  the  Harkov 
approximation  also  Involves  the  neglect  of  terms  of  fourth  and  higher  order 
In  V. 

The  steady  state  (or  equilibrium)  for  the  subsystem  requires  special 

consideration.  That  such  a  state  actually  exists  Is,  for  example,  discussed 
1 4 

by  Van  Hove.  It  Is  characterized  by  J(t)  •  8(«)  being  independent  of 
time.  It  follows  immediately  that  the  Markov  approximation  can  be  made 
without  any  error.  Noting  that  —  -  0  for  any  strength  of  the 
interaction  then  leads  to  the  conclusion  that  the  contribution  to  Eq.  (28) 
of  the  terms  for  each  power  of  V  should  be  zero  independently.  Therefore, 
the  steady-state  solution  can  correctly  be  obtained  from  the  expression 
within  the  Born  and  Markov  approximations  [Eq.  (3D],  even  if  the  Born  and 
Harkov  approximations  are  not  valid  otherwise.  Note  that  then  the 
relaxation  constants  [Eq.  (32)]  do  not  correctly  describe  the  time  evolution 
to  equilibrium,  however  close  the  initial  state  and  the  steady  state  might 
be. 

5.  Comparison  with  reservoir  theory 

An  alternative  procedure  for  obtaining  an  approximate  expression  for 
e(t)  is  reservoir  theory.1*1*  The  starting  point  Is  the  equation  of  motion 
for  the  complete  density  operator  [eq.  (4)]  in  the  interaction  picture, 

^dt^  *  (33) 

where  V(t)  and  P(t)  are  defined  analogously  to  Eqs.  (13)  and  (27). 
Substitution  of  the  formal  solution  of  Eq.  (33)  on  the  right-hand  side  and  a 


change  of  variables  In  the  integral  gives 

t 

-  <j£)tm).p(o)]  ♦  [t(t),[t(t-f>.m-f]].  on) 

o 

An  expression  for  the  tine  derivative  of  B(t)  can  be  obtained  by  taking 
on  both  sides  in  Eq.  (3k),  although  it  still  depends  on  the  coaplete  density 
operator  J(t). 

Equation  (3k)  is  exact,  and  comparison  of  this  equation  with  the 
expression  obtained  with  the  Zwanzlg  projection  technique  [Eq.  (11)]  shows 
that  both  Involve  a  term  depending  on  the  initial  time  density  operator.  In 
reservoir  theory  there  is  a  term  of  apparently  second  order  in  the 
Interaction,  and  dependent  on  the  complete  density  operator.  This  term 
corresponds  to  the  third  term  in  Eq.  (11).  the  latter  which  Includes 
contributions  of  all  orders  in  V,  but  Involves  only  the  subsystem  reduced 
density  operator.  Apparently  the  contributions  of  higher  order  in  V  are 
contained  in  the  coherences  between  the  subsystem  and  the  reservoir,  still 
present  in  Eq.  (34).  The  first  approximation  to  be  made  is  therefore  the 
factorization  of  Mt), 

P(t)  -  pbB(t)  ♦  AJ(t),  (35) 

tdiere  again  the  conditions  of  Eqs.  (9).  (10)  and  (15)  are  imposed  on  pfe  and 

V.  Substituting  Eq.  (35)  into  (3k)  and  discarding  all  terms  involving 

Aft(t),  we  obtain  expression  which  is  identical  to  Eq.  (17).  That  is,  the 

factorization  of  B(t)  in  reservoir  theory  corresponds  to  the  Born 

approximation  in  the  Zwanzlg  projection  technique.  Indeed,  it  can  be  shown 

1 2 

that  the  conditions  for  validity  for  both  approximations  are  the  same. 

Then,  the  Markov  approximation  (replacement  of  B(t-t')  by  8(t))  is  also  the 
aaa»  in  both  approaches.  As  in  the  Zwanzlg  projection  technique,  the 


factorisation  of  9(0)  In  tha  first  torn  [Eq.  (3*)]  introduces  an  error  of 
order  V2,  but  a  factor  ~  saaller  than  the  second  tern.15,1*  The 

h 

factorisation  of  9(t-t')  in  the  second  tern  introduces  an  error  of  order  V  . 

Therefore,  the  neglect  of  Ap(0)  is  an  essentially  different  approxiaation 

than  the  factorisation  of  9(t-t'). 

As  shown  in  Sections  2-A,  the  tern  Involving  p(0)  and  the  error 

introduced  by  the  extension  of  the  Integral  froo  t  to  infinity  are  of  the 

sane  order  in  V  as  the  relaxation  constant,  but  are  saaller  by  a  factor 

t  /t.  Therefore,  the  equation  of  notion  for  9(t),  derived  within  the  Born 
c 

and  Markov  approx i nations  [Eq.  (31)],  la  only  valid  for  tines  t  »  i  . 

c 

Since  the  condition  for  validity  of  the  Born  and  Markov  approxi nations 
[Eqs.  (26)  and  (30)]  requires  the  relaxation  tine  of  the  subsysten  to  be 
■uch  larger  than  the  correlation  tine  of  the  reservoir,  this  tine  Interval 
still  Includes  aost  of  the  transient  regine  before  the  subsysten  reaches  the 
steady  state. 

Disregarding,  for  the  aonent,  the  natrlx  character  of  the  equation  of 
notion  for  9(t)  [Eq. (31)],  the  solution  for  B(t)  is  (tc  «  t  «  T  ') 

T- 

9(t)  -  *Or t  ♦  n  0(-|))B(0),  (36) 

where  the  terns  of  0(-£)  represent  the  neglected  contributions  of  9(0)  and 

the  extension  of  the  Integral.  The  prefactor  of  this  tern  is  essentially  of 

the  sane  order  of  aagnltude  as  the  relaxation  constant  Itself.  Equation 

(36)  shows  that  the  relative  error  in  the  value  of  9(t)  is  of  the  order 

t^/t.  However,  the  relative  error  in  the  different  B(t*t  )  -  9(t)  is  only 
o  c 

2 

of  the  order  (t  /t)  .  Since  the  value  of  this  difference  is  of  the  order 
0 

Te/t,  it  follows  that  Eq.  (31 )  describes  the  variation  of  9(t)  on  any  tine 
scale  with  the  sane  relative  acouracy  as  the  absolute  value  of  9(t)  is 


known.  Therefor*,  the  so-called  coarse-grained  averaging  is  not 

6.  Phonon  relaxation  of  an  adsorbed  molecule 


Thermal  desorption  and  laser -Induced  desorption  have  been  the  subjects 
5—8  17-21 

of  many  studies.  *  The  nodel  employed  by  many  authors  for  the 

vlbrationally  bounded  molecule  is  that  of  a  one-di  mens  Iona 1 

c  9  20  21 

oscillator.  *  '  *  Only  motions  of  the  molecule  perpendicular  to  the 

surface  are  taken  into  account.  This  approach  assumes  that  the  effects  of 

(frustrated)  rotations  of  the  molecule  on  the  substrate  can  be  neglected. 

Although  we  cannot  give  a  rigorous  proof  for  the  correctness  of  this 

assumption,  there  are  some  experiments  which  support  it.  For  example,  for 

CO  on  copper,  it  was  shown  that  CO  adsorbs  in  an  on-top  conformation,  with  a 

relatively  small  angle  (-  10  deg.)  for  the  amplitude  of  the  rotational 
22 

vibration  mode. 

The  interaction  between  the  admolecule  and  the  substrate  is  usually 
described  by  a  Horse  potential,  tftilch  gives  a  fair  description  of  the 
adsorbate-substrate  potential.  It  also  has  the  feature  that  the  eigenvalues 
and  eigenstates  can  be  obtained  analytically,  thus  slmplyfylr^  the  further 
analysis.  The  discrete  levels  are  Identified  with  bound  states  of  the 
adsorbate,  and  the  continuum  levels  correspond  to  gaseous  molecules. 
Desorption  is  described  by  a  transition  from  a  bound  state  to  a  continuum 
state.5 

One  interesting  problem,  which  is  also  experimentally  accessible,  is 
the  calculation  of  the  llnewldth  of  the  vibrational  adbond  due  to  the 
Interaction  with  the  lattice  vibrations  of  the  substrate.  V*  denote  the 


Horse  potential  by  vB(t-tQ),  with  i  the  distance  of  the  admolecule  above  the 
surface  when  the  lattice  vibrations  are  absent.  The  proper  definition  of 


th«  interaction  Hamiltonian  between  the  adbond  vibration  (aubsystea)  and  the 
lattice  vibratlona  (reservoir),  in  accordance  with  Eq.  (10),  is* 

V(*‘VV  "  %(**Vus)  '  <V«o“M«,>'  (37> 

where  <...>  •  Tr  (p®...)  denotes  the  average  over  the  phonons,  and  u  - 
p  o  z 

us(t)  is  the  z -component  of  the  temporal  vibration  amplitude  of  the  surface 

atoa  due  to  the  lattice  vibrations.  The  subsysten  Haelltonlan  is 

«e  -  T  ♦  <VB(*-*0-u2)>.  (38) 

where  T  denotes  the  kinetic  energy  of  the  admoleeule.  The  averaged 

potential  <V(z-z  -u  )>  Is  again  a  Morse  potential,  but  with  renorssllzed 
o  z 

values  for  Its  depth  and  equilibrium  distance,  zQ.*  By  using  Eq.  (37),  the 

relaxation  constants  for  the  adbond  reduced  density  operator  can  be  obtained 

within  the  Born  and  Markov  approximations  in  a  straightforward  calculation.* 

In  the  resulting  expression,  the  reservoir  is  represented  by  the  occurrence 

of  only  the  displacement  autocorrelation  function  <u  (t)u  (o)>. 

z  z 

To  test  the  applicability  of  the  Born  and  Markov  approximations,  we  can 
use  Eqs.  (26),  (29),  and  (30),  tftilch  give  as  the  condition  for  validity 

Yt  «  1.  (39) 

P 

The  characteristic  time  of  the  displacement  autocorrelation  function, 
can  be  identified  with  half  the  phonon  lifetime.  For  metals,  experimental 
determinations  yield  quantities  of  the  order  10  12s. 23,2,1  An  experimental 
determination  of  the  llnewidth  (.Y)  is  available  for  the  0*1  transition  of 
the  B1...C  vibrational  mode  of  CO  adsorbed  on  a  Ni(100)  surface.  A  value  of 
150m"1. 2. 8  1012  a*1  is  found.25  Then 

Yt  -  1. A,  (AO) 

P 

whioh  is  of  the  order  one.  It  follows  that  for  this  system  the  Born  and 
Markov  approximation  are  questionable,  at  best. 


The  estimate  mentioned  above  is  obtained  by  considering  phonons  only 

pC 

(l.a.,  the  phonon  life  tlae  Is  used).  It  was  shown  by  Persson  and  Persson 

that  for  Metallic  substrates  a  da aping  Involving  a  coupling  of  the 

vibrational  aode  with  the  conduction  electrons  will  be  such  larger  than  the 

phonon  damping.  The  strength  of  the  coupling  between  the  vibrational  adbond 

■ode  and  the  electronic  degrees  of  freedoa  Is.  amongst  others,  dependent  on 

the  derivative  of  the  adpartlcle  M.O.  energy  with  respect  to  the  vibration 

coordinate.  For  the  weak  ad -bonds,  as  are  considered  here,  this  derivative 

will  be  much  smaller  than  for  the  Internal  vibration  considered  by  Persson 

and  Persson.  Also,  the  phonon  damping  will  be  much  stronger  for  the  low- 

frequency  adbond  than  for  the  internal  mode.  Therefore,  we  believe  that, 

relative  to  the  electron-loss  mechanism,  the  phonon -damping  will  have  a  much 

larger  contribution  to  the  linewidth  than  in  the  case  of  Internal  vibrations 

of  a  molecule  on  a  metallic  substrate. 

Theoretical  calculations  of  the  relaxation  constants,  performed  within 

the  Born  and  Harkov  approximations,  have  been  published  by  Hood  et  al.27  and 

Volokitin  et  al.  For  CO  absorbed  on  Cu,  a  system  similar  to  CO  on  Nl,  a 

value  of  3.1  lO1^  s  1  was  obtained  for  the  0+1  transition  at  T  •  300  K. 27 

In  this  calculation  a  Debye  spectrum  was  used  for  the  phonon  dispersion 

relation.  For  CO  adsorbed  on  Nl,  a  value  of  Y  •  1.*  lO1^  a*1  was  obtained, 

also  using  a  Debye  spectrum  and  at  T  •  300  K.2®  With  t  •  5  10  1  ^  a,  we 

P 

obtain  Ytp  >  1,  and  It  follows  that  the  Born  and  Markov  approximations  are 

not  Justified  within  the  model  systems  employed  by  those  authors.  The 

values  for  Y  quoted  above  are  obtained  for  an  Infinite  phonon  lifetime.  The 

effect  of  a  finite  phonon  lifetime  Is  to  enlarge  the  calculated  Y  by  at 

27  *1 2 

least  a  factor  ten.  (A  value  t  -  10  s  was  used  in  those 


calculations.)  The  Born  and  Harkov  approx last Ions  are  than  even  more 
questionable. 

Using  a  aore  realistic  surface  phonon  spectrum  due  to  Black  ,  It  was 

possible  to  reproduce  the  experimental  value  for  the  N1...C  vibration  very 

well  (Y  -  13.3  ca  ).  However,  an  Infinite  phonon  llfetlae  was  used. 

Since  a  finite  llfetlae  effectively  broadens  the  phonon  spectrua,  and  since 

29 

the  spectrua  used  Is  sharply  peaked  ,  It  is  to  be  expected  that  inclusion 
of  a  finite  llfetlae  will  give  an  even  more  dramatic  Increase  of  Y  than  for 
the  Debye  spectrua.  Again,  a  much  larger  value  than  experlaent  will  be 
obtained.  The  conclusion  Is  that  the  Born  and  Markov  approx laat Ions  are  not 
Justified  for  this  system. 

CO  adsorbed  on  Nl  or  Cu  is  an  example  of  a  strongly -bound  physlsorbed 

system.  For  the  weakly-bound  physlsorbed  system  Ar  on  W,  a  value  of  Y  - 

11  27 

6x10  s  was  obtained  at  T  •  30  K  and  using  a  Debye  spectrua.  (Note  that 

T  •  30  K  for  Ar  and  T  •  300  K  for  CO  both  correspond  to  the  situation  k-T  < 

»  * 

where  Is  the  fundamental  vibration  frequency.)  Now  Yip  •  0.3,  and 
the  Born  and  Markov  approximation  alght  be  valid,  especially  since  a  aore 
realistic  surface  phonon  dispersion  relation  Is  expected  to  reduce  the 
calculated  value  of  Y  considerably.  Because  the  fundamental  vibration 
frequency  Is  smaller  than  the  Debye  frequency,  the  inclusion  of  a  finite 
In  the  calculation  of  Y  hardly  affects  the  result. 27 


7.  Hu  It 1 phonon  processes 

The  Interaction  potential  between  the  adaolecule  and  the  substrate  can 
be  expanded  In  a  power  series  In  u|t 


(z-x  )  i  »  V*’*n)  2 

V  (x-x  -u  )  •  V  (z-x  ) - r- — — U  ♦  r  - - — ~  --U 

aoz  BO  n  t  2  Jjl  * 


(*1) 


15 


iaa&^i^aaaaiaaiBagaaiMt 


which  shows  that  ths  Interaction  Nasiltonian  can  Indead  be  written  as  a  sub 


of  products  of  s  subsystea  and  s  reservoir  operator.  The  latter  are 

identified  with  the  various  powers  of  u(.  The  subsystea  operators  are  the 

derivatives  of  the  adbond  potential.  va<*~*0)  *•  included  in  the  subsystea 

Haailtonlan.  The  tera  proportional  to  u”  gives  rise  to  n-phonon  processes. 

l.e.,  a  transition  between  levels  of  the  adbond  is  aceoapanled  by  the 

ealsslon  and/or  absorption  of  n  phonons.  It  can  be  shown  that,  for  a  Morse 

potential,  the  contribution  of  n-phonon  processes  to  the  relaxation  constant 

is  exactly  given  by  teras  Involving  <uf(t)uz(0»n  in  the  expression  for  Y 

obtained  by  using  the  full  potential  of  Eq.  (37).^°  For  a  systea  in  which 

the  fundaaental  vibration  frequency  is  larger  than  the  Debye  frequency 

(e.g.,  CO  on  Nl  or  Cu),  the  aultlphonon  processes  give  an  laportant 

contribution  to  the  relaxation  constant.®**7*^0 

It  has  been  argued  that  %*>en  aultlphonon  processes  in  second-order 

perturbation  theory  give  a  significant  contribution,  higher  orders  in  the 

6  7  S 

interaction  also  need  to  be  included.  *  That  this  is  not  the  case 
already  follows  froa  the  fact  that  the  conditions  for  validity  of  the  Born 
and  Markov  approx  last  ions  [Eqs.  (26)  and  (30)]  are  obtained  without  any 
assuaptlon  for  the  fora  of  the  Interaction  potential.  Explicitly,  the 
aultlphonon  contributions  to  Y  are  given  by  [Eq.  (2*)] 


(*2) 


in  second-order  perturbation  theory.  In  fourth-order  of  the  Interaction 
they  are  [Eq.  (23)3 


(«3) 


where  vR  defines  the  contribution  of  n-phonon  processes  to  the  relaxation 


constant.  For  n>1.  2  it  la  given  by 


v 


▼ 


1 


2 


2  »V*-*0>  2 

V*  (V(0)S(-gH  >  1 


a  >v  , 

V.  <vl0)u*(— Jr-11 


(*!**) 

(*»5) 


with  slallar  expressions  for  the  higher  values  of  n.  The  relative 

laportance  of  each  of  the  aultlphonon  processes  follows  from  the  ratio  of 

the  corresponding  i.e.,  it  Is  determined  by  v2/v2.  However,  the 

contribution  of  higher  orders  In  perturbation  theory  Is  determined  by  the 

2  2  2 

smallness  of  the  parameters  v  t  /h  .  More  precise,  they  can  always  be 

n  c 


neglected  If  [Eq.  (26)] 


(Ivn)Tc 
n  n  c 


«  1. 


«6) 


where  the  summation  runs  over  all  the  aultlphonon  contributions.  The 
conclusion  is  that  the  Incorporation  of  aultlphonon  processes  within  second- 
order  perturbation  Is  a  valid  procedure,  only  as  long  as  the  Born  and  Harkov 
approximations  are  Justified  [Eqs.  (26),  (30)  and  (*6)]. 

The  fact  that  Gortel  et  a  1 . ^  find  an  appreciable  contribution  of  the 
fourth-order  terms  aerely  Indicates  that  the  Born  and  Markov  approximation 
cannot  be  aade.  It  then  is  questionable  that  teras  of  up  to  fourth  order  in 
the  Interaction  are  sufficient  to  calculate  Y  accurately.  Moreover,  the 
fourth-order  correction  as  obtained  by  Gortel  et  al.®’7'®,  or  as  given  by 
CO*  C 03 ) .  only  includes  the  teras  connected  with  the  Born  approxiastion. 
Additional  correction  teras  of  fourth  order  In  V  arise  due  to  the  Markov 
approxiastion  not  being  valid  anyaore. 

6.  Coherent  excitation  of  an  adsorbed  aolecule 


An  adsorbed  molecule  Irradiated  by  laser  light  is  a  widely  studied 
system. 17  Examples  of  processes  studied  are  laser-induced  surface  chemical 
reactions  and  laser-induced  desorption.  lhe  effect  of  the  laser  is  to 
directly  excite  the  vibrational  adbond  mode  or  an  Internal  mode  of  the 
adspecles.  Subsequently,  relaxation  occurs  through  the  interaction  with  the 
lattice  vibrations,  as  is  described  in  previous  sections. 

The  system  can  now  be  considered  as  a  subsystem  (the  adbond)  In  contact 
with  two  reservoirs  (the  phonons  and  the  laser,  respectively).  Within  the 
Born  and  Markov  approximations,  a  master  equation  is  obtained,  with 
relaxation  constants  given  as  the  sum  of  a  phonon  part  and  a  term  due  to  the 
laser -ad molecule  interaction, 

y  •  Y  ♦  Y  .  (U7) 

P  r 

The  expression  for  Yp  is  identical  to  the  one  obtained  for  the  adbond  with 

the  phonons  only  [Eq.  (32)].  The  laser  part,  Yp,  can  be  obtained  In  an 

analogous  calculation,  employing  the  laser-adbond  Interaction  as 
20 

perturbation. 

By  an  analysis  analogously  to  Section  3.  It  can  be  shown  that  the 
condition  for  validity  of  the  Born  approximation  Is  the  same  as  If  each 
reservoir  Is  considered  Independently,  we  obtain. 


Vp  <<  1 ' 

<*8a) 

Vr  «  1  * 

(A8b) 

However,  the  condition  for  validity  of  the  Markov  approximation  Is 
different.  It  is  given  by  [Eq.  (29)]: 

Y  t.  «  1 ,  (*9a) 

p  r 


(h  9b) 


That  la.  the  change  of  B(t)  la  required  to  be  saall  on  a  tlaeecele  of  the 
longest  correlation  tine,  tp.  It  follows  that  the  conditions  for  validity 
of  the  Born  and  Markov  approx last Ions  are  not  the  saae,  as  was  seen  to  be 
the  ease  when  only  one  reservoir  la  present  (Sections  3  and  4).  Generally, 
they  are  aore  restrictive  for  the  Markov  approxlaatlon  than  for  the  Born 
approximation. 

Dependent  on  the  relative  strength  of  the  laser  adbond  and  phonon 

adbond  Interaction,  the  conditions  Eqs.  (48)  and  (49)  can  be  simplified.  A 

number  of  Halting  cases  are  given  In  Table  1.  In  Section  6  It  was  shown 

that  for  actual  systems  Y  t  Is  aore  or  less  of  the  order  one.  Since  t  >> 

P  P  r 

t,  it  follows  that  the  conditions  for  validity  [Eqs.  (48)  and  (49)  and 
P 

Table  1]  can  not  be  matched,  except  in  the  case  of  a  weak  laser: 

Y  <  (-2)y.  .  Table  1  shows  that  the  condition  for  validity  of  the  Markov 
r  tp  o 

approxlaatlon  is  the  aost  restrictive  one.  It  is  only  in  this  case,  with 

Y  t  «  1  but  Y  t  «M  ,  that  the  generalized  master  equation  derived  by 

P  P  P  r 

Berl  et  al.10  and  Pereaans  et  al.11  applies.  However,  the  nuaerlcal 

exaaples  given  by  these  authors  refer  to  laser  intensities  stronger  than 

allowed  by  Table  1.  It  Is  interesting  to  note  that  the  Markov  approxlaatlon 

does  not  become  better  when  the  laser  power  Is  reduced  further.  The 

relative  error  In  Y  remains  of  the  order  Yt„ .  This  Is  easily  understood  If 

p  r 

we  realise  that  the  Markov  approxlaatlon  requires  8(t)  to  have  only  a  saall 

ehange  on  a  tlae  r.  This  change  Is  deteralned  by  both  the  phonons  and  the 

r 

lamer  [Eq.  (4?)].  if  y  «  y  ,  then  the  laser  has  relatively  little  effect 

r  p 

on  the  changes  Induced  in  B(t). 

A  situation  of  great  experlaental  Interest  is  when  we  start  with  an 
adbond  in  equlllbrlua  with  the  phonons,  and  then  turn  on  the  laser.  This 
la,  for  exaaple,  encountered  in  spectroscopy .  The  condition  under  which  the 


Born  approximation  la  Justified  is  again  given  by  Eq.  (48).  However,  the 
condition  for  validity  of  the  Harkov  approximation  is  different  now. 
Initially,  there  is  no  change  in  8(t)  due  to  the  phonons.  When  the  laser  is 
turned  on,  an  upper  limit  for  the  change  in  8(t)  over  a  time  Interval  t  is 
given  by  Yr  t.  The  actual  change  will  be  smaller  because  the  phonon 
interaction  forces  B(t)  into  the  direction  of  equilibrium.  The  condition 
for  validity  of  the  Markov  approximation  then  becomes 

Y  t  «  1,  (50) 

r  r  * 

T 

instead  of  Eq.  (49).  It  follows  that  for  weak  lasers,  Y  <  (“*)Y  ,  the 

r  t  d 
r 

condition  for  validity  of  the  Born  and  Harkov  approximations  is  the  same  as 
if  each  reservoir  could  be  considered  independently. 

Finally,  we  want  to  discuss  the  situation  where  the  Born  and  Markov 
approximations  are  made  with  respect  to  one  reservoir  only,  e.g.,  the 
phonons.  For  laser  intensities  which  lead  to  a  change  in  B(t)  smaller  than 
the  change  Induced  by  the  phonons,  the  conditions  under  which  the  Born  and 
Harkov  approximations  are  Justified  are  the  same  as  in  the  absence  of  the 
laser  [Eq.  (30)].  However,  if  the  rate  of  change  of  B(t)  due  to  the  laser 
is  larger  than  the  due  to  the  phonons,  the  Markov  approximation  is  only 
Justified  if  the  change  in  8(t)  in  a  time  due  to  the  laser  is 


9.  Cone 1 us ions 

The  conditions  under  which  the  Born  and  Harkov  approximations  can  be 

■ede  within  the  Zwanzig  projection  foraalisn  have  been  analyzed.  It  is 

shown  that  the  condition  for  validity  of  both  approx  last  Iona  la  the  saae.  A 

ooaparlson  with  reservoir  theory  is  Bade,  and  it  la  pointed  out  that  the 

Born  approx laation  in  the  Zwanzig  projection  technique  corresponds  to  the 

factorization  of  the  density  operator  in  reservoir  theory. 

Apart  from  the  Born  and  Markov  approxiaatlons.  the  derivation  of  a 

master  equation  [Eq.  (3D]  relies  on  two  other  approxiaatlons:  extension  of 

the  upper  bound  of  the  Integral  over  tine  froa  finite  t  to  infinity,  and  the 

factorization  of  the  Initial  tlae  density  operator.  Both  approxi net  ions 

introduce  errors  In  the  relaxation  constant  of  the  order  t  /t,  but  which 

c 

are  of  the  saae  order  in  the  interaction.  Therefore,  the  approxlaate 
equations  apply  only  for  tlaes  t  »  t  . 

C 

For  t  «  T  ,  the  value  of  B(t)  is  proportional  to  V  t,  with  an  error 

due  to  the  initial  tine  factorization  of  the  order  y  t  .  On  a  tlae  scale 

c 

t  ,  9 ( t)  changes  by  an  aaount  V  t  :  however,  the  error  changes  only  by  an 

C  C 

t 

aaount  of  the  order  (r*)Yt  «  Yt  .  Therefore,  the  details  of  the  tlae 

v  C  C 

evolution  of  B(t)  are  preserved  on  any  tlae  scale,  although  its  absolute 
value  is  only  known  with  an  error  of  the  order  Y  tq.  This  conclusion 
contradicts  the  idea  of  a  coarse-grained  average,1’12'15’16  and  indeed  in 
the  derivation  presented  here  a  coarse-grained  averaging  procedure  was  not 

necessary. 

The  results  of  Sections  3*5  are  applied  to  the  case  of  phonon 
relaxation  of  the  vibrational  levels  of  an  adatoa.  It  is  (bund  that  for  CO 
adsorbed  on  HI  or  Cu  (an  example  of  a  strongly-bound  physlsorbed  system), 
the  Born  and  Harkov  approximation  are  not  Justified.  For  V  on  Cu  (an 


21 


example  of  a  weakly-bound  physlsorbed  system),  It  is  shown  that  the  Born  and 
Harkov  approximation  can  be  made. 

In  Sections  6-8  the  Interaction  of  laser  radiation  with  a  phonon-damped 
vibrational  adbond  Is  discussed.  This  system  can  be  seen  as  a  small 
subsystem  (the  adbond)  In  contact  with  two  reservoirs.  It  Is  shown  that  the 
conditions  for  validity  of  the  Born  and  Markov  approximations  are  different 
from  each  other,  and  In  fact,  for  the  Markov  approximation,  more  restrictive 
than  In  the  case  with  a  single  reservoir  [Eqs.  (48)  and  (49)].  For  a  weak 
laser,  l.e.,  one  with  much  less  exciting  power  than  the  relaxation  power  of 
the  phonons,  the  condition  for  validity  of  the  Born  approximation  is 
(Table  1) 

Vb  «  ’■  (51) 

and  the  condition  for  validity  of  the  Markov  approximation  is 

Tkt  «  1.  (52) 

b  r 

It  follows  that  there  Is  a  region  of  strength  of  phonon  interaction  where 
the  Born  approximation  can  be  made,  but  the  Markov  approximation  not.  Then, 
a  generalized  master  equation  should  be  used  to  describe  the  time  evolution 
of  B(t). 10,11  However,  for  an  initial  condition  where  the  adbond  Is  in 
equilibrium  with  the  phonons,  the  equation  for  validity  of  the  Markov 
approximation  [Eq.  (52)]  reduces  to  l.e.,  the  same  as  Eq.  (51).  It 

follows  that  for  analysis  of  spectroscopic  data,  the  Born  and  Markov 
approximations  are  Justified  If  the  reservoir -subsystem  interaction  obeys 
Eq.  (51).  However,  when  spectroscopy  Is  used  in  non-equilibrium  situations 
(m.t* »  to  follow  the  course  of  a  chemical  reaction),  the  reservoir-subsystem 
Interaction  has  to  obey  the  much  more  severe  condition  of  Eq.  (52).  For 
phonon  relaxation  this  can  hardly  be  expected  to  be  the  case. 
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V 


Table  1.  Conditions  for  validity  of  the  Born  and  Markov  approx  last  ions  for 
various  ratios  of  the  strength  of  the  laser  and  phonon 
interactions. 

Batio  of  relaxation  constants  Born  approxlaatlon  Harkov  approximation 
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